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Abstract 

To investigate an effect of social interaction on the bystanders' intervention in 
emergency situations we introduce a rescue model which includes the effects of the 
victim's acquaintance with bystanders and those among bystanders. This model 
reproduces the surprising experimental result that the helping rate tends to decrease 
although the number of bystanders k increases. The model also shows that given 
the coupling effect among bystanders, for a certain range of small k the helping rate 
increases according to k and that coupling effect plays both positive and negative 
roles in emergencies. Finally we find a broad range of coupling strength to maximize 
the helping rate. 
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1 Introduction 



Recently the concepts and methods of statistical physics and nonlinear dynam- 
ics are applied to investigate social, economic and psychological phenomena 
[1,2,3]. Among the interesting subjects that have attracted physicists are the 
opinion dynamics [4,5] including voting process [6,7] and social impact the- 
ory [8,9,10,11,12,13]. Social impact theory stemmed from the bystander effect 
that was first demonstrated in the laboratory by psychologists who studied the 
effect of social interaction among bystanders [14] . This effect is a marvel phe- 
nomenon where people are less likely to intervene in emergencies when others 
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are present than when they are alone. A well known example is Kitty Gen- 
ovese case [15]. She was stabbed to death at one night in 1964 by a mentally 
ill murderer for about 30 minutes while at least 38 witnesses were believed 
to have failed to help the victim. On the other hand the refutation that the 
popular account of the murder is mostly wrong was suggested on the web- 
site [16]. Regardless to this controversy [17] many experimental results in the 
laboratory and from field studies show the bystander effect obviously [14,18]. 

Although the social impact theory originated from the bystander effect and 
related models have been developed, there does not exist any model describing 
quantitatively the bystander effect itself. To construct the model mentioned 
above we abstract a few important factors from the existing literatures [14,19]. 
The main dependent variables are the rates of reporting the smoke in the test 
room, helping endangered victims and so on. When tested alone, subjects 
behaved reasonably and the response rate was high. However the rate was sig- 
nificantly depressed when they were with other subjects. Subjects with others 
were unsure of what had happened or thought other people would or could do 
something. In another experiment subjects who were friends responded faster 
than those who were strangers. The subjects who had met the victim were 
significantly faster to report victim's distress than other subjects. And the 
degree of arousal that bystanders perceive is a monotonic positive function of 
the perceived severity and clarity of the emergency, and bystander's emotional 
involvement with the victim, which will be also considered by our rescue model 
in an abstract way. 



2 Rescue Model 



A simple rescue model is introduced for investigating the statistical feature 
of the effects of the victim's acquaintance with bystanders and those among 
bystanders in emergency situations. We focus on the relations between agents 
rather than agents themselves, so define a relation spin between two agents as 
which is 1 if agents i and j have been successful in the intervention and 
otherwise, aij can be interpreted as an element of adjacency matrix of helping 
network. Each agent i has its intervention threshold q over which that agent 
can try to intervene in an emergency situation. 

At each time step an accident happens which consists of the degree of accident 
represented as a random number q v uniformly drawn from [0, 1], a randomly 
selected victim v, the number of bystanders k confined to [1, N — 1], and a 
set of randomly selected k bystanders N v from a population. Then the update 
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rule is as following: 



a vi (t+l) = e lq v + aa vi (t)+f3 ^ (2a„(t) - 1) - cA (1) 

where 9{xi) is a heaviside step function. Only one bystander i e N v with the 
largest value intervene per accident. If we assume that the response 

speed of bystander i is exponential in x^, the selection of the bystander with 
the largest Xt is justified. Additionally, once one bystander intervenes, the 
pressures on the others will disappear. If Xj > 0, the rescue succeeds and 
then the victim v and that bystander i gain acquaintances if they have not 
been related to each other. In case of Xi < the rescue fails and then their 
acquaintance is cut if existed, a represents the degree of victim's acquaintance 
with bystander, so can be called an acquaintance strength. The third term of 
9 function is related to the acquaintances among bystanders, for each relation 
2a,ij — 1 gives 1 if two bystanders know each other or —1 otherwise. There does 
not exist any neutral relation here. (3 is used to tune the strength of coupling 
so can be called a coupling strength. As an order parameter we adopt the 
average helping rate: 

= N(N-l) S^ift) ( 2 ) 

which can be also interpreted as a social temperature or an average linkage 
of the helping network. We fix q = c = 0.25 for all i according to the ex- 
perimental result [14] that 70 ~ 75% of alone subjects intervened and c does 
not change through this paper, which means we consider a population com- 
posed of homogeneous and non-adaptive agents. And for the most cases a is 
also fixed as 0.1. The main control parameter is the number of bystanders k. 
Finally, the initial conditions are = for all pairs. 



2. 1 Case with (3 = 



At first let us consider only the effect of victim's acquaintance with bystanders 
and ignore those among bystanders. Generally an equation for the average 
helping rate can be written as 

^ = WW " W^ (3) 

where dk(t) denote (a)(t) for a fixed k. Given the values of parameters as 
above, when < q v < c — a, the rescue fails independent of a V i(t) so one of 
acquaintances, if exists, should be cut. When c — a < q v < c, a vi (t+l) = a vi {t) 
so this interval of q v does not contribute to the equation. When c < q v < 1 the 
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rescue succeeds independent of a V i(t) so a new acquaintance is formed only if 
there is no acquaintance between them. 



W ->i = (l-c)(l-a k (t)) k (4) 
W^ = (c-a)(l-(l-a k (t)) k ) (5) 

The stationary solution is obtained easily. 

"* = 1 -(t^) < 6 > 

The solution means that although the coupling strength is not considered 
((3 = 0), the helping rate depends on the number of bystanders k. As increasing 
k the probability that the victim does not know any bystanders (contributes 
to Wq^i) decreases rapidly and the complementary probability (contributes to 
Wi^q) increases. This is why the helping rate decays. In this case a affects not 
the decaying behavior but the decaying speed only when a < c. If a becomes 
larger than c, Wi^o = so a k = 1 independent of k. 

The numerical simulations are shown in Fig. 1 and reproduce the experimental 
results that people who live in larger cities are less likely to help the strangers 
[18]. Compared to Ref. [18] k corresponds to the size of city or town. 



2.2 Case with (3 > and k = 2 



If the coupling among bystanders is taken into account and k = 2, the interval 
[0, 1] of q v is divided into five subintervals by four points c — a — (3, c — a + (3, 
c — (3, and c + (3. Whether c — a + (3 is smaller or larger than c — (3 does not 
affect the result. For each subinterval the transition rate is calculated. 



W ^i = (1 - c - (3)(l - a 2 f + (1 - c + (3)a 2 (l - a 2 ) 2 (7) 
W^ = (c-a + (3)a 2 (l - a 2 )(2 - a 2 ) + (c - a - (3)a\(2 - a 2 ) (8) 

The second term in R.H.S. of W\^q vanishes when c — a — (3 < 0, so we get 
the solutions for both cases. 



a 2 



l-a-/3- v /a 2 +/3 2 +a/3-a-/3+c(l-a) 

1-Q 

l+c-2a-/3- v /5c 2 +4a 2 -3/3 2 -8ac-2/3c-2c+2/3+l 



\ic-a- (3>0 

(9) 



2(c-a-/3) 



ii c - a- ft <0 



The numerical simulations in Fig. 2 support the analytic solution. One can 
find a maximum value of a 2 when c — a — (3 = 0, which means that there exists 
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Fig. 1. The analytic solution and numerical simulations of the case with a > and 
P = 0. Here N = 100, c = 0.25 and a = 0.1. 

an optimal coupling strength (3 opt = c — a to maximize the helping rate when 
given c and a. 

2.3 General case with f3 > 

For a general k, the interval [0, 1] of q v is divided by the following points: 
c — a — (k — 1)0, c-a+(k- 1)(3, c-(k- 1)(3, • • • , c+ (k— 1)13. For 
each subinterval the transition rate can be calculated. 

Wo->i = (1 - a k ) k (1 - c + (k - i)/3 - 2(3F(a k , k)) (10) 
W^ = (l - (1 - a k ) k ) (c - a - (k - l)(3 + 2(3F{a k , k)) (11) 

where 

F(a k , k) = |> - 1 - n) f ^ (A; " 1} <(1 - a,)^^ 1 )-" (12) 

n=0 \ n I 

Consequently to get a stationary solution we must solve the following equation. 

(1 - a)(l - a k ) k — c — a — (k — l)[3 + 2f3F(a k , k) (13) 
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Fig. 2. The analytic solution and numerical simulations of the case with a, (3 > 
and k = 2. a2 has a maximum value when c — a — (3 = 0. The other parameter 
values are the same as Fig. 1. 

Since F(a,k, k) is the |fc(fc — l)th order polynomial, the order of above equation 
is max(|A;(£; — 1), k) and difficult to solve exactly for all k. Therefore, as an 
approximation we consider only the most left, the second and the most right 
subintervals. 



Wo-,i = (l — c—(k — 1)8) (1 - a k ) k (14) 

W^ = ( c -a-(k- 3)8) (l - (1 - a k ) k ) (15) 

The approximate stationary solution is 

/ c — a — (k — 3)8 \ , c — a . . 

a * =l - ( l- a -2(t-2W " i " t£hS T U <16) 



Figure 3 shows the numerical simulations supporting the approximation ex- 
cept for k Pi k\. The helping rate undergoes a trough for k less than fci, 
which is also partly observed in the experiment results [18], and then de- 
creases monotonically as for the case without a coupling effect (8 = 0). For 
k > A; 2 the helping rate shows a sudden drop to 0. To understand this be- 
havior intuitively two factors should be considered. On the one hand as for 
the case without a coupling effect, the increasing number of bystanders in- 
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Fig. 3. The analytic solution and numerical simulations of the case with (3 = 0.01. 
Circles represent the numerical results and solid line and dotted line do the approx- 
imate solution for f3 = 0.01 and for f3 = for comparison respectively. Inset graph 
shows the exponentially decaying behavior of the same data. 



creases the probability of cutting acquaintance too. On the other hand the 
most left subinterval [0,c— a — (k — 1)13), which most contributes to W\^q, 
decreases as increasing k. In this paper c is fixed as a rather small value so 
that a change of the most right subinterval rarely affects 111 conclusion, 

when k < k the former factor dominates the latter so that the helping rate 
decreases. When k Q < k < k\ the increasing k minimizes the inhibiting effect 
by other bystanders to make the helping rate increase. Once k passes ki, the 
decreasing subinterval vanishes and only the former factor works. 



Coupling effect among bystanders plays both positive and negative roles in 
emergencies. When k < k , the helping rate is rarely affected by the coupling. 
In case of ko < k < k 2 , the coupling helps the helping rate obviously. Especially 
when k\ < k < k 2 , decays exponentially with slope of —0.0268 in the 
semi-log plot. Once k passes k 2 , turns to because of too many other 
inhibiting bystanders. The value of k 2 can be obtained by a simple mean field 
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Fig. 4. The maximum value a max as a function of (3. There exists a broad range of 
(3 to maximize the helping rate. 



approximation and by using the following: 



(6(x - x )) x 



1 if x < 
1 — x if < x < 1 
if x > 1 



(17) 



means an average over a; G [0, 1]. From Equation (1), if c — aa k — f3(k 
l)(2afc — 1) > 1, then = 0, therefore the condition for = is k > k 2 



— + 1 



Under what conditions is the helping rate maximized? For each (a, /?) we can 
find the maximum value of a k , denoted by a max . It is sufficient to compare a\ 
to a kl , where k\ is defined by Eq. 16. a-y is ^— ^ independent of (3 while a kl 
is smaller than a\ for small /5 and exceeds a\ for some range of (3 and then 
vanishes, in fact k± — 1 for larger /? in Fig. 4. Especially a maa; approaches to 1 
when 0.08 < (3 < 0.35. There exists a broad range of (3 for which the helping 
network is almost fully connected. 
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3 Conclusions 



We introduced a simple rescue model to investigate the effects of victim's 
acquaintance with bystanders (acquaintance strength a) and those among 
bystanders (coupling strength 0) for the bystanders' intervention in emergency 
situations. When (3 = 0, as the increasing number of bystanders k the helping 
rate decreases, where the speed depends on a. For the case of (3 > and k = 2, 
there exists an optimal coupling strength (3 = c — a for the maximum helping 
rate. Coupling strength plays both positive and negative roles in emergencies. 
For ko < k < k\, since the coupling among bystanders minimizes the inhibiting 
effect, the helping rate increases according to k. And then the helping rate 
decays monotonically as for the case without a coupling effect. Once k passes 
k 2 too many bystanders inhibit the helping. There exists a rather broad range 
of f3 where almost all the trials to intervene in emergencies are successful. 
In this case it is not necessary to fine-tune the coupling strength to get the 
highest helping rate. 
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